arXiv:1509.01946v2 [math.DG] 25 Feb 2016 


Implicit Lagrange-Routh Equations and Dirac Reduction 


Eduardo Garcia-Torano Andres* 

Department of Mathematics, Faculty of Science 
The University of Ostrava 
30. dubna 22, 701 03 Ostrava, Czech Republic 

Tom Mestdagl 
Department of Mathematics 
Ghent University 

Krijgslaan 281, B-9000 Gent, Belgium 
Hiroaki Yoshimura^ 

Department of Applied Mechanics and Aerospace Engineering 

Waseda University 
Okubo, Shinjuku, Tokyo 
169-8555, Japan 

26 February 2016 

Abstract 

In this paper, we make a generalization of Routh’s reduction method for Lagrangian 
systems with symmetry to the case where not any regularity condition is imposed on the 
Lagrangian. First, we show how implicit Lagrange-Routh equations can be obtained 
from the Hamilton-Pontryagin principle, by making use of an anholonomic frame, and 
how these equations can be reduced. To do this, we keep the momentum constraint 
implicit throughout and we make use of a Routhian function defined on a certain sub¬ 
manifold of the Pontryagin bundle. Then, we show how the reduced implicit Lagrange- 
Routh equations can be described in the context of dynamical systems associated to 
Dirac structures, in which we fully utilize a symmetry reduction procedure for implicit 
Hamiltonian systems with symmetry. 
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1 Introduction 

There is no doubt that there exists a close relation between symmetries and conservation 
laws, which has been one of the fundamental motivations for many geometric approaches 
to mechanical systems. The symmetry group of a dynamical system can always be used to 
reduce the system to one with fewer variables. When symmetry, besides, leads to conserved 
quantities, it can be very advantageous to incorporate that property into the reduction 
process. For example, when the system is Hamiltonian on a symplectic manifold, one first 
restricts the attention to the submanifold determined by the conserved momenta, and only 
later one takes the quotient of this submanifold by the remaining symmetry (which in general 
happens to be only a subgroup of the original symmetry group). This, in a few words, is the 
so-called symplectic reduction theorem (see, Marsden and Weinstein [1974]; Marsden [1992]). 

Symplectic reduction may be applied to the standard case of classical Hamiltonian sys¬ 
tems defined on the cotangent bundle. While this procedure has been thoroughly studied in 
the literature, its Lagrangian counterpart, the so-called Routh or tangent bundle reduction, 
has traditionally received much less attention, even though since its conception in Routh 
[1877, 1884] it has proven to be a valuable tool to obtain and discuss, e.g., the stability of 
steady motions or relative equilibria. A few modern approaches to the topic can be found 
in the papers Crampin and Mestdag [2008]; Langerock, Cantrijn and Vankerschaver [2010]; 
Marsden, Ratiu and Scheurle [2000]. One of the drawbacks of these papers is that a regu¬ 
larity condition needs to be assumed. In this paper, we will focus upon Routh reduction 
within the context of Dirac structures, without assuming any regularity hypotheses. 

For simplicity, let us consider for a moment the case of a Lagrangian L{x,±,6) with a 
single cyclic coordinate 9. The first step in Routh’s procedure is to write the corresponding 
velocity 9 in terms of the remaining coordinates and velocities (x, x) by making use of the 
conservation law dL/d9 = /r, which follows from Noether’s theorem. One then introduces 
the restriction R^{x, x) of the function L — 9{dL/d9) to the level set where the momentum 
is the so-called Routhian (see, e.g., Marsden [1992]). With this function, one can observe 
that the remaining Euler-Lagrange equations of the coordinates x, again when constrained 
to the level set associated to are in fact Euler-Lagrange equations for the Routhian R^^. 
The end result of Routh’s reduction method is therefore that it reduces the Euler-Lagrange 
equations of the Lagrangian L(x, x, 9) to those of i?^(x, x) on a reduced configuration space. 
A crucial ingredient in the above process, however, is that the Lagrangian satishes the reg¬ 
ularity condition {d'^L/d^9) ^ 0, which is necessary for carrying out the first step. Routh’s 
procedure and the regularity condition for it to be applicable can be generalized to arbi¬ 
trary Lagrangians with a (possibly) non-Abelian symmetry group G. In this situation, the 
condition is often referred to as G-regularity. 

It is easy to construct a Lagrangian which fails to be G-regular. The following example 
in is taken from Langerock and Castrillon [2010]: 

L(x, y, Vx, Vy) = {vxf' -f VxVy -V{x). 

Note that it has a cyclic coordinate y (and therefore an Abelian symmetry group G = R), but 
that it is not G-regular. Also linear G-invariant Lagrangians will always fail to be G-regular. 
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For example, the dynamics of N vortices in the plane admit the following Lagrangian: 

^k^k') 2 EE 7 „ 7 fcln|z„ - Zfe|, ziGC, 

k n k^n 

where 7 fe £ K are parameters of the model; see Chapman [1978] for more details. This 
Lagrangian is clearly linear in its velocities and invariant under rotations of the vortices in 
the plane but not G-regular. Also in the context of plasma physics, linear Lagrangians often 
appear (see, e.g., Littlejohn [1983]). 

The aim of this paper is to extend Routh’s method to the most general case where 
not any regularity condition is imposed on the Lagrangian. Our approach is based on the 
Hamilton-Pontryagin principle (as it is called in Yoshimura and Marsden [2006b]) which 
leads to an implicit formulation of the Euler-Lagrange equations on the so-called Pontrya- 
gin bundle TQ®T*Q. We will show that under the assumption of symmetry, we can reduce 
these implicit equations to a set of reduced implicit Lagrange-Routh equations. The key 
ingredient is that we can circumvent the hypotheses on regularity by keeping the momen¬ 
tum constraint implicit throughout. Our method involves a generalized Routhian function 
which is defined on a certain submanifold of the Pontryagin bundle rather than on a sub¬ 
manifold of the tangent bundle, as is commonly the case for G-regular systems. Implicit 
Lagrangian systems can be geometrically described in the framework of Dirac structures 
(see Yoshimura and Marsden [2006a]). The definition of Dirac structure in Courant [1990]; 
Dorfman [1993] was originally inspired by the notion of Dirac brackets, which was coined 
by Paul Dirac (in the 1950s) for dealing with constraints in the Hamiltonian setting when the 
given Lagrangian is singular (see e.g. Courant and Weinstein [1988]; van der Schaft and Jeltsema 
[2014]). So from the very start, there has been a strong relation with singular Lagrangians 
and constraints. In the second part of the paper, we will show that the reduced implicit 
Lagrange-Routh equations may be also formulated in terms of a Dirac structure, by con¬ 
sidering a reduction method known for implicit Hamiltonian systems from van der Schaft 
[1998]; Blankenstein and van der Schaft [2001]. 

For completeness, we mention that the paper by Langerock and Castrillon [2010] also 
deals with the general case. However, these authors use a variational approach which is 
based on Hamilton’s principle rather than on the Hamilton-Pontryagin principle. Therefore, 
it focusses on different aspects of the theory. 

This paper is organized as follows. In §2, we review the derivation of the standard 
implicit Euler-Lagrange equations for a possibly degenerate (or singular) Lagrangian L via 
the Hamilton-Pontryagin principle. We use a technique that is similar to the one that has 
been used in, for instance, Crampin and Mestdag [2008, 2010] to rewrite the implicit Euler- 
Lagrange equations in terms of an anholonomic frame. Once the implicit equations on a 
general frame are obtained, we specialize these expressions to a particular frame adapted 
to a given symmetry of the Lagrangian (§3). For a prescribed value of momenta, we find 
the implicit Lagrange-Routh equations, and express them in an invariant form. In §4 we 
reduce them to obtain the reduced implicit Lagrange-Routh equations. The regular cases 
are discussed in §5, where we illustrate our theory by showing how the reduced implicit 
Lagrange-Routh equations agree with those developed in the literature. §6 rephrases the 
previous results in terms of reduction of Dirac structures. We show how the reduced implicit 
Lagrange-Routh equations correspond to a certain reduced implicit Hamiltonian system. 
Finally, in §7, some examples are shown. 
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2 A version of the Hamilton-Pontryagin principle using 
anholonomic frames 


Hamilton-Pontryagin principles. Let Q be a configuration manifold of a mechanical 
system with dimQ = n. Coordinates on Q are given by q°‘, fiber coordinates on TQ 
and T*Q will be denoted by v°‘ and pa, respectively. In the following, the index a runs 
from 1 to n unless otherwise noted. The notations are chosen in such a way that we can 
make a notational difference between a general curve {q{t),v{t)) in TQ, and the lifted curve 
{q{t), q{t)) in TQ of a curve q{t) in Q, where t denotes the time in I = {t £ R \ a < t < b}. 

The Hamilton-Pontryagin principle leads to an implicit form of the Euler-Lagrange equa¬ 
tions of a possibly degenerate Lagrangian L. These equations follow from considering the 
following variational principle on TQ 0 T*Q: 

[L{q,v) + {p,{q-v))]dt = S [ [L{q°^,v°‘) + Pc{q°'- v°‘)]dt = 0, 

Jb Jb 


for variations of {q{t),v{t),p(t)) where q{t) has fixed endpoints and v{t) and p{t) are ar¬ 
bitrary. From this, we can easily conclude that a solution {q{t),v{t),p(t)) of the implicit 
Euler-Lagrange equations must satisfy 




n 


Pa = 


dL 


( 2 . 1 ) 


See Yoshimura and Marsden [2006b] for more details. 


Anholonomic frames and quasi-velocities. In this section, we shall rewrite the im¬ 
plicit Euler-Lagrange equations in terms of the so-called quasi-velocities. Lagrangian equa¬ 
tions which involve quasi-velocities are often called Hamel equations in the literature (see, for 
instance, Marsden and Scheurle [1993] ; Bloch, Marsden and Zenkov [2009] ; Crampin and Mestdag 
[2010]). We will need these expressions when we consider the Routhian in the following sec¬ 
tions. 

In the next paragraphs, we will need the natural lifts of vector fields on Q to its tangent 
manifold and Pontryagin bundle, respectively. Let {q°‘,v°‘) be the natural tangent bundle 
coordinates on TQ. li X = X°'{d/dq°‘) is a vector field on Q, then its complete lift X^ and 
vertical lift X^ are the vector fields on TQ, given by 


X® = ^ + 


d dX^ 


d 


dqP dq'T dvP' 


x^ = xP 


dv^ 


Likewise, its complete lift to M = TQ 0 T*Q is the vector field 


= Y" 


dq<^ 


dXP d dXP d 


A standard reference for the properties of these vector fields is the book Yano and Ishihara 
[1973]. Given a vector field Y = Y°‘{d/dq°') on Q, we can form the linear function Y = Y'^pa 
on T*Q C M. Likewise, for a 1-form 0 = O^dq'^ we can define a linear function on 0 = 9aV°‘ 
on TQ CL M. The following properties can then easily be verified: 


Y^(Y) = [Y, Y], X^{9) = Lx9. 


( 2 . 2 ) 


Quasi-velocities are fiber coordinates in T^Q, defined with respect to a non-coordinate 
or anholonomic frame. Let Za = Z^{d/dq^) be a new basis for the set of vector fields 
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on Q. This means that at each point q the matrix {Z'^(q)) has an inverse matrix, smoothly 
defined, which we will denote by {Wp{q)). Each vector Vq € TqQ can be expressed by 
Vq = v°‘Za{q). The fiber coordinates (v“) are then the quasi-velocities of Vq with respect 
to the frame {Za}. Their relation to the natural fiber coordinates is simply v“ = . 

The coordinate frame {d/dq°‘} is an example of a frame, whose corresponding quasi¬ 
velocities are simply the natural fiber coordinates v°‘. A measure for the deviation of a 
given frame {Za} from being a coordinate frame, is given by its object of anholonomity 
(see, e.g., Schouten [1954]), which is defined by the relation 


[ZfjjZy] — R’p^Za- 

The are given in coordinates by the following expressions: 


DO! _ 

- 


dZi 





= -{z 




dq^ 


(2.3) 


We can lift the frame {Za} on Q to the frame {2’®, Z'^} on TQ. In what follows, we 
will often make use of the following, easily verifiable, properties: 

zC(/) = zf, ^I(/) = 0, 

^a(v^) = -<v^ = (2.4) 

The 1-forms tT“ = W^dq^ form a basis for 1-forms on Q. If {q°‘,Pa) denote the natural 
coordinates on T*Q, we can also introduce quasi-momenta by means of p /3 = Z^pa, and 
we note that the natural pairing is preserved: {p,v) = PaV°^ = Pav“. 

The independent variations Jpa} form a basis for all variations on M = TQ © 

T*Q. We can change this to a new basis, adjusted to the new frame {Za^ IE“} on M. If we 
denote 

= W^q^, 

then {r(;“, 5v“, (5 pq} will be all independent variations and hence this set will form a new 
basis for all variations on M. Here one should think of the quasi-velocities v“ as functions 
on TQ (or on M), and therefore 

dq'y ^ 


The implicit Lagrangian systems with quasi-velocities and quasi-momenta. The 

direct computation using (2.3) yields the variation of a Lagrangian L on TQ as 

or or 

SL = —6q^ + —Sv^ = ZS{L)w^ + Z^{L){6v^ + (2.5) 

We will now give a version of the Hamilton-Pontryagin principle that makes use of the 
anholonomic frame. We are looking for a curve (^“(t), v“(t), Pa(t)) in M = TQ © T*Q, 
namely the one whose base curve is q{t) and whose fiber coordinates are given by the curve 
(v“(t), Pq)^)) in quasi-velocities and quasi-momenta which satisfies the following variational 
principle 

pCL pa 

0 = 6 [L{q,v)-i-{p,q-v)]dt = 6 [T(g, v) + Pa(u“ - v“)] dt. 

Jb Jb 

Here u'^(t) stand for the quasi-velocities of the lifted curve q{t) in TQ, namely, u'^(t) = 

W^iqm^t). 
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If we take the above expression (2.5) for SL into account, we obtain 
0 = / [Z^{L)w° + Z^{L){Sv° + + 5p„(u“ - v“) + pa5u“ - pa5v“] dt. 

Jb 

First, we compute 


/ o pa oW 

p„5u“dt = + 


dt 


lb L 


dq'i dq!^ 
[{PaR'/lr'^^ - Pr)w^] dt. 


- PaW 


dt 


In the above, we have used integration by parts and the fact that 6q°‘{a) = 5q°‘{h) = 0. 
We have also made use of the expression (2.3) for R’^^.. Then, by substituting this into our 
variational principle, it follows that 


0 = ^ [(Z‘;{L) + pci?“,u'^ + Zl (L)i?“^v^ - p,) 

+ (^Z^ (L) - Pa) (5v“ - (^V“ - U“) (Spo 


dt 


holds for all variations w'^, 6v°‘ and (5pa. Since these are all independent, we conclude that 
the curve (^“(t), v“(t),Pa(t)) should satisfy 

v“ = u“, Pa=^^(T), Pa = ^°(T). (2.6) 


These are the implicit Euler-Lagrange equations in quasi-velocities and quasi¬ 
momenta. Remark that the last equation takes such a simple form, only because we have 
already made use of the first two equations. 


3 The implicit Lagrange-Routh equations 

The Lie group action on Q. Assume that C? x Q —>• Q is a free and proper action of 
a possibly non-Abelian Lie group G on Q, so that we can regard Q —>■ Q/G as a principal 
G-bundle. We will assume in this section that the Lagrangian A ; TQ —> R is invariant 
under the symmetry group G. We use cc* for coordinates on the shape space Q/G and 
((?“) = {x^,9°‘) for coordinates on Q. We denote the corresponding coordinates on the 
tangent bundle by As before, we use the notation {x{t),x{t)) to denote the lifted 

curve on T{Q/G) of a curve x{t) in Q/G. 

We now choose a principal connection on Q —Q/G and let Xi denote the horizontal lift 
of the coordinate vector fields d/dx^ on the shape space Q/G with respect to this connection. 
In terms of the coordinates introduced above, we have 

X= — -A^ — 

* dx^ * 50 “’ 

where A“ are functions on Q (the connection coefficients) and we employ the symbol A to 
denote the connection one-form on TQ which takes values in g. 
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Let {Ea] be a basis of the Lie algebra g, and {-E®} the corresponding dual basis of g*. 

We denote by the structure constants of g, [Ea,Eb\ = C^j^Ec- The fundamental vector 
fields associated to the action will be denoted by {Ea}- We can express them as 

F = K^ — 

for some functions on Q, often called the coefficients of the infinitesimal generator map. 

We suppose throughout the paper that G is connected. This has the advantage that 
invariance of functions and tensor fields can be checked by the vanishing of the Lie derivatives 
of these functions and tensor fields, in the direction of the fundamental vector fields of the 
action. 

The action on Q lifts to actions on TQ, T*Q and M = TQ (B T*Q. In each case, it is 
well-known that the fundamental vector fields of the lifted actions are given by the complete 
lifts (to TQ, T*Q and M, respectively) of the fundamental vector fields on Q. In the case 
of the action on M the fundamental vector fields are therefore (linear combinations of) the 
vector fields E^. Since G is supposed to be connected, a function E on M is invariant if, 
and only if, E^ (E) = 0. 

The implicit Lagrange-Routh equations. Let us rewrite the implicit Lagrange equa¬ 
tions (2.6) by making use of a specific anholonomic frame. First we consider the frame 
{Za} = {Xi, Ea} on Q. This corresponds to the so-called moving frame in literature. We 
denote the corresponding quasi-velocities and quasi-momenta by (v\v“) and (pi,pa). In 
fact, since the vector fields Xi are assumed to project onto the coordinate fields on Q/G, 
the quasi-velocities v* can be naturally identified with the natural fiber coordinates v’’ on 
T{Q/G). The curve u*(t) that appears in the equations (2.6) is then simply the lifted curve 
i*(t) of the curve x'^{t) in QjG. Similarly, the quasi-momenta pi can be identified with the 
momenta pi of T*{Q/G), but this identification is not canonical since it depends on the 
choice of the connection A. From now on, for simplicity, we will use the notation u® and pi 
to denote the corresponding quasimomenta. 

The brackets of the frame are given by: 

[X„ X,] = E“ Ea, [X„ Ea] = 0, [Ea, Efc] = -C:,Ea, 

where the E“- stand, up to a sign, for the curvature coefficients of the principal connection A 

(this is the convention in Crampin and Mestdag [2008], but differs from e.g. Marsden, Ratiu and Scheurle 

[2000]). The relation \Xi,Ea] is a consequence of the invariance of the Xi. 

The Lagrangian is invariant if, and only if, E^ (L) = E^{L) = 0. The implicit Euler- 
Lagrange equations (2.6) are therefore 

v“ = u“, Pa=Ej(E), Pa = 0, 

E=±\ p, = X7(L), p,=XpiL). 

From the top row, we see that Pa is constant along solutions, say f>a = Fa, with p = 

PaE‘^ G g*. We can define a generalized Routhian, as the function on TQ given by 

R^{q,v) = L{q,v) - pa^°‘- ( 3 . 1 ) 

It is, however, not the standard definition of the Routhian function as one may find in, for in¬ 
stance, Marsden, Ratiu and Scheurle [2000]; Crampin and Mestdag [2008]; Langerock, Cantrijn and Vankerschaver 
[2010]. We will clarify the relation between these two definitions later in in §5. The advan¬ 
tage of the current definition is that it allows us to keep the momentum constraint implicit. 
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From the relations (2.4) of the previous section, we obtain 

Xf{R>^) = Xf{L) + (L), 

Ej{R^ = {L) - 

Therefore, a solution of the implicit Euler-Lagrange equations (2.6) is a curve 
(x*(t),r(t),u*(t),v“(t),p,(t),pa(t)) ■.IcR^XI = TQ®T*Q 


satisfying 

E^{Rn = 0, Pa=^ia, 

XY{Rn = Pm P.= X^{Rn - PaBt.vE 

We will call these equations the implicit Lagrange-Routh equations. The terminol¬ 
ogy Lagrange-Routh equations is adopted from Marsden, Ratiu and Scheurle [2000] . 

The implicit Lagrange-Routh equations in invariant form. We will restrict our 
attention to one specific level set of momentum. Consider the canonical momentum map 
J : T*Q —?► g* of the G-action on Q, and fix a value g G g* . Let XI ^ denote the submanifold 
TQ 0 J~^{g) in M = TQ 0 T*Q. Local coordinates on are then {q,v\v^,pi), the Pa 
being fixed by the value g G g* . 

The G-action on M restricts to a G^-action on M^, where G^ stands for the isotropy 
group. We will describe how solutions of the implicit Lagrange-Routh equations (3.3) which 
happen to lie on can be projected to curves in M^/G^^, satisfying some reduced equa¬ 
tions. In order to do that, we need to rewrite them in such a way that all involved terms 
are given by G^-invariant functions. When that is the case, these G^-invariant equations 
will project to equations on M^/G^. 

Let 5 = ^°'Ea € g be an arbitrary element of g, then it follows from (3.2) that = 

—^“^cG^^v^. From this we see that ,^®(i?^) = 0, if and only if, ^ G g^. Therefore, we see 
that R^ is (only) G^-invariant. The Routhian R^ can thus be identified with a reduced 
function on TQ/Gfj, (for which we shall use the same notation). 

We next define local coordinates on M^/G^. It is easy to see that the quasi-velocities u® 
and the quasi-momenta pi are G^-invariant functions on M^. Indeed, from (2.2) and (2.4), 
we have _ 

Ef(«®)=0, Gf(p,) = [4,^.] = 0. 

The last property is based on the observation that Xi = pi. The above expressions show 
that and pi (thought of as coordinate functions on M) are G-invariant, and thus also Gf^- 
invariant functions on M (and therefore also on M^). On the other hand, the quasi-velocities 
v“ are not G^-invariant functions on since 

(v^) = -Glr. 

To overcome this issue we introduce a new frame that is completely G-invariant (and 
therefore also G^i-invariant), see also Crampin and Mestdag [2008]. This coincides in the 
literature with the so-called body-fixed frame. Consider a new set of vector fields, given 
by Ea = The following reasoning shows that there exists a matrix {A^) of functions 

on Q for which these vector fields are all invariant. They may be invariant if and only if 
0 = [Ea,Eb] = (^EaiAl) — C^^A'^'^Ec- The integrability condition that is needed for the 
PDE equation 

EfiAl) - Glb^Ai = 0 


v“ = u“, 
u® = x\ 
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to have a solution A'^ is satisfied by virtue of the Jacobi identity of the Lie bracket on 
0 . We can therefore claim that, at least locally, the above PDE has a solution for which 
A = {A\) is non-singular, and for which A is the identity on some specified local section of 
TT : Q ^ Q/G. 

An explicit way to define the vector fields Ea, and the one we will use henceforth, 
is as follows. Let U C Q/G be an open set over which Q is locally trivial. Then the 
fibration is n : U x G ^ U, and the action is given by ijjg{x,h) = {x,gh). We can define 

Ea ■■ {x,g) AdgE^{x,g) = Tipg{Ea{x, e)), where Adg : 0 —0 is the adjoint action. In 
coordinates, we write 

F = T^ — 

where L\ are functions on Q. 

Another way to think of these two frames is the following: If Q is the Lie group G, and 
the action is given by left multiplication then the tilde-vector fields coincide with a basis of 
right-invariant vector fields, while the hat-vector fields are all left-invariant. 

The quasi-velocities (v*,v“) with respect to {Xi,Ea} are all invariant functions on M, 
since now 

Af(v*')= 0 . 

We can therefore take , z;®, v“,pi) for our local coordinates on M^/G^, where [q]G^ 

stands for the coordinates of the orbit of q under the G^-action. In a more global inter¬ 
pretation, we have split the quotient Mg^/Gg, by making use of the principal connection A 
as 

Mg/Gg ~ {Q/Gg) Xq/Q (T (Q/G) © 0 © ( Q/G) ) , 

where § is the adjoint bundle. In fact, the coordinates above correspond to fiber coordinates 
with respect to this identification (see Langerock and Castrillon [2010] for more details). 

We now check whether all the terms that appear in the implicit Lagrange-Routh equa¬ 
tions (3.3) are G^-invariant. For example, the function E^{R^) is not G^-invariant, but 
the function E^(R^) = (R^) is- Indeed, for ^ G 0 ^, 

^E^ [E^m) = {eE^m) = 0 , 

where we have used that [E^,E^] = [Eb,Ec]‘^ = 0 , and where we have also used that the 
Routhian is G^-invariant. We conclude that we need to replace the equation E^{R^) = 0 
in (3.3) by the equivalent equation E^(R^) = ^b^a (R^) = 0 to obtain a G^-invariant 
equation. With a similar argument we can show that the functions Xp{R^) and XY{R^) 
are all G/^-invariant. 

Remark that, since [Xi, Xj] is an invariant vector field, we must have that [A^, BfjEY] = 
0 , or equivalently 

EYiBY) - BtgCl = 0. 

Using this, we can now check that the function g,aBfjV^ is G^-invariant; namely, for ^ G 0 ^, 
we find 

eEYifiaBYv^) = = 0 . 

To conclude, apart from the defining relation pa = g,a for Mg, the following system of 
equations is equivalent to the implicit Lagrange-Routh equations (3.3) and consists only of 
G;_j-invariant equations: 

r=n“, EY{Rn = 0, 

E=x\ xYiRn=P^, 


p, = xY{Rn-FaBYvE 


(3.4) 
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All these equations correspond to equations on whose coordinate expressions 

will be obtained in the next section. 

We finish this paragraph with a coordinate expression for the equations v“ = u“. Recall 
first that u“ is the quasi-velocity corresponding to the lifted curve {q,q). If we write 

we can conclude that the equation v“ = u“ is equivalent with the equation v^L^ = 


4 Reduction of the implicit Lagrange-Routh equations 


We shall use the residual G^-symmetry of the equations (3.4) to drop them to M^/G^. 

To do this, we will make use of an invariant decomposition (with respect to the adjoint 
action of G^) of the Lie algebra fl = meaning that, for all g £ we have 

Adg(fl/fl;, ) C (g/O/i). In terms of coordinates, this splitting will be denoted as follows: we 
choose a basis {Ea] = {Ea^ A/} of g in such a way that {Ea} represents a basis of g^ and 
{Ei} a basis of (g/g^x)- 

An invariant splitting g = g^x © (g/g^x) of the Lie algebra as above, together with a 
principal connection on the bundle Q —>■ Q jG, induces a principal connection on the bundle 
Q QjG^. Indeed, let denote the projector TTg^ : g —g^x with respect to the previous 
decomposition. Define the principal connection 1-form A^ := o A : TQ —s- g^, where 
A : Q —?> g is the principal connection 1-form on the bundle Q —?> QjG. This is a well-defined 
connection 1-form because 

(1) If C e g^, then certainly = Tx^^{^)=t,\ 

(2) If g € G^, then since Adg is linear and Adg(g^) = g^, it follows that TTg^ o Adg = 

Adg OTTg^. 

Therefore, A^ is G^-equivariant. Although we will not explicitly make use of it, we men¬ 
tion that an (Ehresmann) connection on Q/G^x —> Q/G can be directly obtained by project¬ 
ing the connection A (note that it is well-defined due to the equivariance of A). This connec¬ 
tion, together with A and A^, plays a role when looking at variational principles in the con¬ 
text of Routh reduction (see also Langerock and Castrillon [2010] and Marsden, Ratiu and Scheurle 
[ 2000 ]). 

Using that g^ is a Lie subalgebra, we have G^^ = 0. On the other hand, from the 
definition of g^, we also get G^i^fic = 0. Finally the fact that the term (g/g^) is invariant 
leads to the relation G^ = 0. We will need these relations later in the paper. 

We will use coordinates (0“) = (9^,6^) such that the fibers of G —G/G^ are given by 
9^ = constant. Then, there are functions on Q such that 


Ea = K^ 


d 


El = Kf 


d 


+ Kj 


d 

W' 


We remark that, by construction, AT^ = 0. Likewise, we can set 


Ea = L^ 


d 

d9B^ 


Ei = Lf 


d 

d9B 


+ Lj 


d 

89-^' 


Since Xj is a G;x-invariant vector field on TQ, it can be thought of as a vector field 
on TQ/G^. This reduced vector field will be completely determined by its action on the 
coordinate functions {x^,9^, u®, v“), which define the coordinates on TQ/G^. We have, again 
using (2.2) and (2.4), 


a7(x®) = 0, x7(0^)=O, XJ{E) = 5], A7(v“)=0. 
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The reduced vector field of Xj is therefore the vector field djdv^ on Likewise, the 

reduced vector fields of are d/dY°‘. It also follows that 


= 5], 


Xf{0^) = -Al Xf{v^) = 0, Xf (v“) = -B 


jfc' 


where [Xi,Xj] = BfjEa (by construction, we have BfjL\ = or equivalently Bfj = 

A’^B\^). Therefore the reduction of X'^ to a vector held on TQ/G^ is 

dxi ^ del Jk Q^a- 


Summarizing, we have proved the following: 

Proposition 4.1. A curve {x''{t),9^{t),v^{t),v°-{t),pi{t)) in M^^/G^ is a solution of the 

reduced implicit Lagrange-Routh equations if it satisfies: 


x‘ = u*, e^ = Lj — x^Aj, 

_ dRf^ dRf^ _ ^ 

~ dv^ ’ 9v“ “ 


^ dRt^ ^jdRf^ 

~ ~ 'W 


e-aBi.x^, 


(4.1) 


None of the above equations depends explicitly on 9^. Therefore, the above equations 
determine the reduced curve on M^/G^. 

Once we have solved these reduced equations for (a;*(t), 9^ {t), v^{t), v“(t),pi(t)) on M^/G^, 
we can recover the complete solution {x^{t),9^{t),9^{t),v^{t),v°-{t),pi{t),pa{t)) on M by 
solving the reconstruction equations 


9^ = - x^Af, Pa = ga 


The hrst equation can only be given a geometrically concise interpretation if we make use 
of a principal connection on the principal bundle TQ ^ TQ/G^, but we will not go into 
the details of this procedure (we refer the interested reader to Crampin and Mestdag [2008] 
for a similar situation in the case of standard Routh reduction). 


5 Special cases 

In this section, we shall obtain the implicit Lagrange-Routh equations in some particular 
cases and relate the resultant equations, when possible, to the results derived elsewhere in 
the literature on Routh reduction. 

The regular case. Let us consider the case where the Lagrangian is regular with respect 
to the group variables. More precisely, a Lagrangian is said to be G-regular if the Hessian 
{Ea {E^{L))) is non-singular everywhere on TQ, which in coordinates can be expressed as 

" {Me:) * 0 ^ 

G-regularity is probably the weakest condition that the Lagrangian should satisfy to have 
an analogy with the classical procedure of Routh. In this case, one of the implicit Euler- 
Lagrange equations (3.4), namely the relation E^{R^) = 0 (or, equivalently, E)^{L) = pa) 
can locally be rewritten in either one of the following explicit forms v“ = or v“ = 

(q, u*), where are smooth functions of (g, z;*). This defines a submanifold of TQ, 







5 Special cases 


12 


with inclusion We can introduce the new function o on This is the func¬ 

tion which is commonly called the Routhian (see, e.g., the papers Crampin and Mestdag 
[2008]; Langerock, Cantrijn and Vankerschaver [2010]; Marsden, Ratiu and Scheurle [2000]). 
It is easy to see that, for its reduced version on N^/G^, we obtain 

d^_(d^ \ (dl^ \ dll 

dx^ \ dx’’ V dv°' dx^ 


In view of the fact that dR^/dv°' = 0 is part of the reduced implicit Lagrange-Routh 
equations (4.1), every instance of dR^/dx'^ in (4.1) can be replaced by dR.^/dx'^, and simi¬ 
larly for dR^ /dv''. The remaining reduced equations are then simply given by 


X = V 


ai 


= li{qy)K^j-x^K 


dt 


^R^^\ 
dv^ ) 


dRt^ jdRf^ 

dx"^ * 89^ 




(5.1) 


The above equations can be found in Crampin and Mestdag [2008]. 

An even stronger regularity condition is the one used in Langerock, Cantrijn and Vankerschaver 
[2010]. Define for each Vq £ TqQ a map S'l ■ fl 0 * as follows: 

3^ - Jl (vq + ■ 

In the above, Jl ■ TQ —^ g* is a standard momentum map on TQ defined by Jl = J o FL, 
where J : T*Q —>■ g* is the standard momentum map. If we assume that is a diffeo- 
morphism for every Vq £ TQ, it has been shown in Langerock, Cantrijn and Vankerschaver 
[ 2010 ] that it is possible to realize the previous equations as the symplectic reduction of 
the original Lagrangian system. Simple mechanical systems, for which the Lagrangian is 
of the form L = T — V with T given by a Riemannian metric on Q, satisfy automat¬ 
ically this stronger form of G-regularity (this follows easily from positive-definiteness of 
the metric). For a detailed study of the Routh reduction for simple mechanical systems, 
see Marsden, Ratiu and Scheurle [2000]. 


The Abelian case. When the group of symmetries is Abelian we have = G. In this 
case, there are no coordinates 9^ and we can just write 0“ everywhere. The equations for 
the curve (x®(t),u®(<), v“(t),pi(t)) become 


^R^^ 


x = V 


Pi = 


p^ = 


dRf^ dRf^ 


= 0 . 


dx^ dv^ ’ (9v“ 

Note that even when the group is Abelian, there remains a curvature term in the equations. 
From the reconstruction equations 


9“ = v''L? - i^A?, 


Pa — Mci? 


we can determine 0“(t) and In the case where L is G-regular, the manifold can 

be identified with T{Q/G) and the reduced equations (5.1) can be regarded as the Euler- 
Lagrange equations of the Routhian on T{Q/G) subjected to a gyroscopic term arising 
from the curvature of A. 

An important case of an Abelian symmetry group is that when the Lagrangian has cyclic 
coordinates. In this case we have a configuration manifold that is a product Q = S x G, 
and the Lagrangian is assumed to be invariant under the action of G on the second factor. 
The equations for (x*(t),r!®(t), v“(t),Pi(t)) become then 


Pi = 


dj^ 

dx^ ’ 


dR>^ dR>^ 


Pi = 


dv^ 


d<7°- 


= 0 . 


Again, if we have G-regularity, these might be further simplified to yield the Euler-Lagrange 
equations of R^ (with no gyroscopic term). 
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6 Routh-Dirac reduction 

In this section, we shall show how the reduced implicit Lagrange-Routh equations can be 
obtained as a reduced Lagrange-Dirac dynamical system. We will call this type of reduction 
Routh-Dirac reduction. The procedure we follow relies on some results known for implicit 
Hamiltonian systems. Again, we assume the same setting as before: a free and proper action 
of a Lie group G on Q, which leaves the Lagrangian L invariant. 

Dirac dynamical system. We refer to Courant [1990]; Dorfman [1993] for the orig¬ 
inal works on the notion of a Dirac structure and to Dalsmo and van der Schaft [1999]; 
Yoshimura and Marsden [2006a]; Cendra, Ratiu and Yoshimura [2015] for more details on 
Dirac structures and the dynamics associated to them. A linear Dirac structure on a vector 
space y is a subspace Dy of H © H* which is a Lagrangian subspace with respect to the 
pairing 

(((r;i, ai), ('C 2 ,0:2))) = (0:2, wi) + («!, V 2 )- 

A Dirac structure on a manifold M is a subbundle Dm CTM (BT*M such that DMim) C 
TjnM X T^M is a linear Dirac structure on T^M at each m G M. In what follows, we will 
use the terminology Dirac dynamical system, as in Cendra, Ratiu and Yoshimura [2015], to 
refer to a wide class of implicit Lagrangian or Hamiltonian systems that can be defined in the 
context of Dirac structures. Given an energy form ip G T{T*M), the dynamics of the Dirac 
dynamical system {ip, Dm) is given by the following condition on a curve c : I cM. ^ M: 

c(t) © p (c(t)) G DM(c(t)). 

Implicit Euler-Lagrange equations. Using the above, we can obtain the implicit Euler- 
Lagrange equations, by taking M to be the Pontryagin bundle M = TQ © T*Q over Q. 
In order to construct Dm we consider the pullback Dm = of the canonical 

symplectic form Qt*q = —d0T*Q on T*Q to M by using the canonical projection ttt*q ■ 
M —>■ T*Q. Then, Dm is a presymplectic form on M, which naturally defines a Dirac 
structure Dm on M by means of the graph of Dm: for each {q,v,p) G M, set 

DMiq,v,p) = {{{dq,Sv,Sp), {a,l3,"f)) G TM ®T*M \ a + Sp = 0, j = Sq, /3 = 0}. 

The generalized energy function EL{q,v,p) = {p,v) — L{q,v) on M defines the energy form 
When we express the two-form Dt*q in natural coordinates as dq°‘ A dp a it is easy to 
check that a curve c{t) = {q{t),v{t),p{t)) in M is a solution of the implicit Euler-Lagrange 
equations (2.1) if (d£i,DM) satisfies the so-called Lagrange-Dirac dynamical system 

c{t) © d£i(c(t)) G DMicft)), for all t. (6.1) 

On the other hand, when we write the generalized energy in quasi-velocities and quasi¬ 
momenta as £l = Pav“ — L and d£i = v“dpQ, + Padv" — dL, and likewise the two-form 
Dt*q as 

Dt.q = 1U“ a dp„ + A (6.2) 

we may easily obtain from (6.1) the implicit Euler-Lagrange equations with quasi-velocities 
and quasi-momenta (2.6). Note that 1U“ has been introduced before as a 1-form on Q, 
while in the above we think of 1U“ as a (semi-basic) 1-form on T*Q. We will use this slight 
abuse of notation from now on. 
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Reduction of Dirac structures. In this paragraph we assume that M is an arbitrary 
manifold, endowed with a Dirac structure Dm- We first recall some generalities on the 
reduced Dirac structure induced by a Lie group action. The results here were originally 
developed in van der Schaft [1998]; Blankenstein and van der Schaft [2001]. 

Assume that a Lie group G acts (freely and properly) on the manifold M, and denote by 
4>g(jn) = g ■ m the action of 5 G G on a point m G M . This action lifts naturally by tangent 
and cotangent lifts to TM 0 T*M. Assume also that the action admits an equivariant 
momentum map relative to Dm-, that is, assume that there exists a G-equivariant map 
: M ^ g* such that |0 dJ^ G Dm, for all C S g, where we recall that | stands for the 
infinitesimal generator of the action associated with ^ G g and is the smooth function 
on M defined by J^{m) = (J^(m),0 , m G M. 

We assume now that D is invariant under G. The Dirac reduction procedure will be 
carried out in the following two steps. First, if /r G g* is a regular value of , then 
C M is a submanifold. If the vector subspace DM{rn) n x 

C TmMfj, X T^M\mij, has constant dimension at each m G M^, then these vector 
spaces naturally induces a restriction of the Dirac structure Dm C TM^ 0 T*M^. 
Second, one observes that the Dirac structure is G^-invariant since 

DM^{g -171)= g- DM^{m), 

where Gf_i = {g € G \ Ad* g = g} is the coadjoint isotropy subgroup of g- This leads to 
a reduced Dirac structure Dm^/g,^ C T{Mfi/G^) 0r*(M^/G^) on the reduced space 
M^/G^ = J~^{g)/Gfj,, which is given by 

Dmjg, ■- = {iX,a) G X{M^/G^) x U}{MJG^) \3{Y,P) G Dm, 

such that Ttt^ oY = X o 7r*a = /3}, 

where tt^ : Mf^/G^ is the canonical projection, which is a surjective submersion. 

Symmetry reduction of implicit Hamiltonian systems. An important case of Dirac 
dynamical systems is that when the energy section is given by the differential of a Hamilto¬ 
nian function H on M. Let H he a G-invariant Hamiltonian on M and let c{t) be a solution 
curve for the implicit Hamiltonian system (diL, Dm), i-e. a curve that satisfies 

c{t) 0 dH{c{t)) G DM^cit))- 


Then, 


-^(c(t)) = {dJ^ {c{t)),c{t)) = - (dH,l^ {c{t)) = 0, for all t and C G 0- 

So, is a first integral of the implicit Hamiltonian system and we can restrict the Dirac 
dynamical system (diJ, Dm) to (diJ^,DM^), where = D|m^- Next, we can reduce 
the restricted implicit Hamiltonian system (diJ^,DM^) to obtain the reduced implicit 
Hamiltonian system (dij^,D m^/g^)i which satisfies 

c{t) 0 dij^(c(t))) G Dm^/gM^), (6-4) 

for each c(f) = 7r^(c(t)) in M^/G^, where -k = is the reduced Hamiltonian on 
M^fG^. 

More details on this reduction of Dirac structures and its associated reduced implicit 
Hamiltonian systems can also be found in Blankenstein and Ratiu [2004]. 
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The reduced implicit Lagrange-Routh equations. We consider again the Dirac struc¬ 
ture Dm on the Pontryagin bundle M = TQ(BT*Q given by the graph of the presymplectic 
form Dm = For : M ^ g*, we take = J o where J : T*Q —>■ g* 

stands for the standard momentum map on Q. The goal of this paragraph is to demonstrate 
that, in this particular setting, the reduced implicit Hamiltonian system (6.4), with Hamil¬ 
tonian H = £l, is nothing but the system given by the reduced implicit Lagrange-Routh 
equations (4.1). 

It is well-known that if Dm is a Dirac structure given by the graph of a symplectic form 
Dm the reduced Dirac structure Dm^/g^ may be given by the graph of a reduced symplectic 
form ^m^/g^- The following observations may be obtained: 

i) The action of G on M = TQ (B T*Q restricts to a G^-action on by tangent and 

cotangent lifts, and this action leaves the presymplectic form Dm^ invariant. Indeed, 
since Dm^ = **Dm (where i : —>■ M is the inclusion), its invariance follows directly 

from the G^-invariance of Dm- 

ii) Moreover, one can show that the form Dm drops to M^/G^. Indeed, it suffices to check 

that Dm^ annihilates vectors which are vertical to the fibration tt^ : —)• M^/G^. 

Thus the reduced presymplectic form ^Mij,ig^ is defined on M^/G^ from the invariance 

of ■ 

iii) It follows immediately that Dm is G^-invariant. In particular, we can define a reduced 
Dirac structure Dm^/g^ on M^/G^, 

Dmjg, c r(M^/G^) ©T*(M^/G^) 


such that 

'^li^M^/G^ = I'*Dm- 

Note that the above characterization of Dm /g can also be obtained in the standard 
cases of symplectic and presymplectic reduction by the well-known characterization of 
the reduced (pre)symplectic form, to be found in Marsden and Weinstein [1974] and 
Echevern'a-Enn'quez, Muhoz-Lecanda and Roman-Roy [1999]). We plan to investigate 
whether the same result may also hold for almost Dirac structures with regular distri¬ 
butions. 

To write down the reduced system of the implicit Hamilton system (6.1), we first need 
to give a coordinate expression of the two 2-forms Dm^ and Dm^/g^ whose graph define the 
Dirac structures Dm^ and Dm^j^/g^- 

Using a principal connection A on the bundle ■. Q ^ QjG, one identifies J~^{pL) ~ 
T*{Q/G) y^Q/G Q- Under this identification the presymplectic form reads 

= ^Q/G - 

where we have used a slight abuse of notations and omitted the pullbacks from the spaces 
T*{Q/G) and Q where the forms ^q/g and are defined, respectively. Erom Cartan’s 
structure equation, it follows dA^ = (^, [A, H]) — Therefore, in coordinates, we have 

Dm^ = dx'‘ A dpi + (^B^jdx'- A dx^ - G^^& A , (6.5) 

where stand for the dual of the basis {Xi,Ea} and where we have identified 

X* = dx^ in the notations of the previous sections. We point out that, again, we do not 
write explicitly the pullbacks: one should think of the forms £1“ in (6.5) as semi-basic forms 
on T*Q or, equivalently, as the vertical lifts {E°‘)^ to T*Q of the corresponding forms 
in Q (see Yano and Ishihara [1973] for more details). To ease the notation we will keep 


6 Routh-Dirac reduction 


16 


this convention from now on since the coordinate expressions agree, unless there is risk of 
confusion. The expression (6.5) is a particular instance of the expression (6.2) in the current 
frame. 

Let be the generalized Routhian given in (3.1) and consider the energy section de¬ 
termined by the differential of the energy E-l restricted to M^. Writing in terms 

of the Routhian as £^ = and using a formula similar to the one we had for 

calculating variations of L in terms of the anholonomic frame (see (2.5)), it follows 

d£^ = - (Xf (i?'^) + dx^ - {XY{Rn - P.) dE 

+ (jia + eJ(R n) - eJ (i?^)dv“ + Edp,, (6.6) 

where we have used the relation E^{R>^) = —pcC^i,v^ from (3.2). 

Recall that C 0 T*M^ is the Dirac structure induced from Before 

we continue with the expression for /q , it is instructive to have a look at the first 
step in the reduction, namely the restriction of the dynamics to M^. Consider for that 
reason the Lagrange-Dirac dynamical system (d£,, ), which satisfies, for each cii) = 
(x*(t),6»“(t),u*(t),v“(t),p,(t)) in M^, 

c{t) © d£^ (c(t)) e ■ 

It leads to the following set of equations: 

X^(R'^)=0, XY{R^^)-P^ = Q, 

PaBt.x^ +p, = Xf (i?^) + E^{RnBt,v\ (pa + {Rn) u^ 

where u^ = (X“^)^(0“ + i;®A“). The above equations might be further simplified as 

PaC^^r = paCS^n^, E^ {Rn = 0, 

E = x\ xY {Rn =K, ^ {Rn - ^^aB:y. 

The similarity with the implicit Lagrange-Routh equations (3.3) is obvious, with the 
only difference that it is not possible to conclude from equations (6.7) that W = u'^. This is 
a consequence of the fact that solutions of a presymplectic equation are only determined up 
to elements in the kernel of the presymplectic form (see Gotay and Nester [1979]). Indeed, 
considering the splitting 0 = © (fl/0/i) introduced earlier in §4, the relation PaC^ci^^ ~ 

u'^) = 0 implies that v'^ = u^, but it is not true in general that also v^ = u"^. This is 
reminiscent of the standard case of symplectic reduction on T*Q; if : J~^{p) —>■ T*Q 
denotes the inclusion, then the kernel of the presymplectic form = z*Dq is given by 
kerD^ = {E \ E & 0^}. 

The equations (6.7) only refer to the restriction step in the reduction process of an 
implicit Hamilton system. As explained before, in the second step, we need to reduce that 
system by G^. This boils down to looking at the Dirac structure defined by the graph of 

Proposition 6.1. Let be the reduced energy function on Af^/G^ defined by = £^. 

Given the Dirac structure a curve c(t) = 7r^(c(t)) in M^/G^ is a solution curve 

of the Routh-Dirac dynamical system (d€^, which satisfies 

c{t) © d£^ (c(t)) e (Dm^/gJ (c(t)), 
if and only if the reduced implicit Lagrange-Routh equations (4-1) hold. 
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Proof. Using the splitting {Ea) = {Ej,Ea) of g, we write the presymplectic form in 
as 

= dx^ A dpi + (^B!^jdx'' A dx^ — CfjE^ A E'^'j . 

Recall that, with the notations in §3, we have Ea = A\Ei,. Then we have, for duals, 
E^ = AjE"^. With this we can rewrite the previous expression of Um in terms of the forms 
E^: 

= dx^ A dp^ + (S“ dx* A dx^ - A^^AiCfjE^ A E^^ . 

In this frame, the reduced two-form formally the same expression. 

The proof will follow from the expression (6.6) for d£^ in invariant (hat) coordinates, 
and from its reduction to This computation is completely analogous to the one carried 
out before in §4. If we take into account that RA = A'^B’^j, = AjV'^, etc., it follows that 


d£, 


fd^ 


V dx^ 


del J 


dx* — 


\ dv^ 


-pAdE-paA^KAiC^adE^E^ 


dvf- 


ri 

C^aAA-^tA^KAi<^^E^ - -^dr + Edp.. 


Using the relations for the structure constants derived in §4, it is then straightforward to 
check that a curve 


c(t) = {x\t),e\t),v''{t),Y°'{t),pi{t)) : / C M ^ 


satisfies (6.4), if and only if the implicit Lagrange-Routh equations (4.1) hold. One of the 
equations will appear rather as paCfj{v^ — u^) = 0, but, as we mentioned before, this 
implies v^ = u^. I 


Remarks. We will call the reduction process that we have just described Routh-Dirac 
reduction. It is a particular instance of a reduced Dirac dynamical system. One can see 
a similar discussion on the Dirac reduction associated with the symplectic reduction in the 
more general context of Dirac anchored vector bundle reduction in Cendra, Ratiu and Yoshimura 
[2015]. 

It should be remarked that the above Dirac structure Dm^/g^ is different from the one 
used in Yoshimura and Marsden [2007, 2009] , where the Dirac structure is rather defined as 
a subbundle of the vector bundle (T M®T*M)/G. This last definition may be advantageous 
when one wants to give a geometric interpretation of the so-called implicit Lagrange-Poincare 
or Hamilton-Poincare equations in the variational link with Dirac structures. These two sets 
of equations are the result of a reduction process that takes the full symmetry group G of a 
mechanical system into account, but it does not produce conserved quantities p. Our aim, 
in the previous sections, was a Routh-type reduction; so we wanted to take full advantage 
of the conserved quantities, at the price of reducing by a possibly smaller symmetry group 

G^. 


7 Examples 

To conclude, we will now discuss two illustrative examples where the regularity conditions 
on the Lagrangian fail. 
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Linear Lagrangians. A linear Lagrangian is probably the easiest case where L is not 
locally G-regular. In general, a linear Lagrangian L on TQ conveys to the form L = 
{a{q),Vq) — f{q) for some 1-form a on Q and function / on Q. We say a few words about 
the two examples we mentioned in the Introduction. 

The Lagrangian 

L{x, y, Vx, Vy) = (vx)^ + VxVy - V(x) 

of the first example has a cyclic coordinate y. Fix a value y, G M. for the momentum py. The 
Routhian, constructed with respect to the trivial connection on R, reads 

R^{x,y,Vx,Vy) = [vxY + VxVy - yvy - V{x), 


and the reduced implicit Routh equations are then 


X=Vx, Px = -^ = -V'{x), 


Px = 2Vx + V, 




— y = 0. 


Together with the reconstructions equations y = Vy and py = y, we get a solution of the 
original implicit Euler-Lagrange equations for L with the prescribed value of the momentum 

Py 

Consider again the Lagrangian for the dynamics of point vortices in the plane from 
Chapman [1978]: 

L{zuk) = ^Y. lk{zkh - ZkZk) - ^ EE 7„7feln|2:„ - Zk\. 

k n k^n 


Although it is invariant under rotations on C, the coordinates are not adapted. A possible 
way to proceed is to take polar coordinates for each position Zk{t) G C, say Zk = pke'^^'", and 
then to consider the relative angle with respect to 6i. Defining (pk = for fc > 2 and 

(pi =01, we have an S'^-action along (pi with associated momentum dL/dpi, where L is the 
Lagrangian in the new coordinates. A computation then shows that this conserved quantity 
is precisely the moment of circulation (also called angular impulse in Newton [2001]) 


I = ^lkpl = ^lk\zk?- 

k k 


From here the implicit Lagrange-Routh equations follow without further difficulty. 


A degenerate Lagrangian. We will now discuss a mechanical model for field theories 
to be found in Capri and Kobayashi [1982, 1987], although we will follow the exposition 
in Echeverrfa-Enriquez, Muhoz-Lecanda and Roman-Roy [1999], where it appears in the 
context of presymplectic reduction. The Lagrangian is 

L = {'ipym^j{ipy + {yycij{tpy - {tpyciyipy - {tpyfiytpy, 

where represent the scalar complex fields (which are regarded as coordinates in the 

model), and where the matrices , Cij , satisfy: rriij , are hermitian and (cij ) = 

—Cji- This guarantees that L is real. 

We will consider the same values for these matrices that appear in Echeverrfa-Enrfquez, Muhoz-Lecanda and Rom 
[1999]: 

/O 0 0 \ ^ /O 0 0\ /l 0 0\ 

m^j = \ 0 7712 0 , dj = c^j = - \0 i 0 , Tij = \0 1 0 , 

\0 0 TO3/ ^ \0 0 ij \0 0 1/ 
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and which lead to a degenerate Lagrangian. With that choice, it is apparent that L becomes 
invariant under the actions of 5"^ by rotation in the second and third scalar fields. Thus, we 
have a T^-action on Q. Writing + iy^ and + iv^ the Lagrangian becomes: 

L = 1712 [{u ) + (v ) j + m3 [{u ) +(r')j+t!a: +vx — u y — u y 

- + iv^f) - + {y^f) - + ( 2 /^)^) ■ 

In this case, coordinates adapted to the T^-action are simply the usual polar coordinates 
on both {x^,y‘^) and We will denote these sets of polar coordinates by (r, 0) and 

(p, p) respectively. Then, it follows 

L = m2 {r'^Vg + vl) + m 3 {p^vl + v^) + r^ve + - r'^ - p^ - {{x^Y + : 

and the equations for pg and are of the form 

pg = 217127^ Vg + r^, P 0 = 2m3p^Vg + p^. 


The Routhian reads, for a given choice {pg-,p^) = {pg, p^), and with the trivial connec¬ 
tion, as follows: 

Ri^ = {r'^vg + vl) + m3 {p^Vg + vl) + r^vg + p^v^, - r'^ - p^ - ((x^)^ -|- (j/^)^) - pgvg - p^v^. 
The reduced implicit Lagrange-Routh equations in this case are 


x^ = u^. 

Px^ 

= —2x^, 

Px^ 

= 0 , 

II 

Pyl 

= -2y\ 

Pyl 

= 0 , 

r = Vr, 

Pr 

= 2rvg + 2rvg, 

Pr 

= 2v,, 

II 

Pp 

= -h 2pv^, 

Pp 

= 2v, 


2r‘^ve+r‘^ = pg, 2p'^v,p + p^ = p^, 

and the reconstruction equations are given by 

0 = vg, P = v^, pg = pe, p 0 = 
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